Magnetodielectric Response of the Spin-Ice Dy 2 Ti 2 07 



Masafumi Saito/'Q Ryuji Higashinaka, 1 and Yoshiteru Maeno 1,2 

'Department of Physics, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan 
2 International Innovation Center, Kyoto University, Kyoto 606-8501, Japan 
(Dated: February 2, 2008) 

We report the magneto-dielectric response of single crystals of the spin-ice compound Dy2Ti2C>7 
down to 0.26 K. The dielectric constant under zero magnetic field exhibits a clear decrease reflecting 
the development of the local two-spins-in, two-spins-out structure below about f .2 K. Both the real 
and imaginary parts of the dielectric constant under magnetic fields sensitively respond to various 
changes in the spin structures. We found that the real part can be described in terms of local spin 
correlations among the moments of tetrahedra, rather than among individual Dy 3+ moments. Using 
the peaks in the imaginary part, we have constructed a precise field-temperature phase diagram in 
the [111] field direction. We thus demonstrate that the magneto-dielectric response can be a high- 
sensitivity local probe of the spin state of geometrically frustrated systems. 



PACS numbers: 75.50.Lk, 75.80.+q, 77.84.Dy 

I. INTRODUCTION 

Geometrically frustrated spin systems are now clearly 
recognized as a class of magnets exhibiting unusual be- 
havior such as residual entropy, spin chirality, and quan- 
tum spin-liquid statei. For the experimental investiga- 
tion of frustrated mechanisms leading to such behavior, 
it is crucial to probe the local changes in the spin config- 
uration not always probed readily by conventional bulk 
methods. 

Recently, the strong coupling between dielectric and 
magnetic properties of titanates^ and manganites^ at- 
tracted renewed research interest. A phenomenologi- 
cal expression for the magneto-dielectric response, which 
takes into account of local spin correlations, has been 
proposed 2 . This motivated us to investigate how the di- 
electric constant responds to a variety of changes in the 
spin configuration in the spin-ice magnet, an archety- 
pal system with geometrical spin frustration. In this 
study, we use the spin-ice compound Dy2Ti2 07 in the 
pyrochlore structure^, in which Dy 3+ ions constitute a 
three-dimensional network of corner-shared tetrahedra, 
so called the pyrochlore lattice shown in Fig.^ Because 
of crystal field splitting of the J = 15/2 multiplets, the 
ground-state of the Dy 3+ ion is a doublet with the mag- 
netic moment exhibiting strong "Ising anisotropy" along 
the local <111> direction, which point to the center of a 
tetrahedron from a vertex. Dipolar interaction dominates 
the nearest-neighbor interaction and is ferromagnetic, al- 
though the antiferromagnetic exchange interaction is also 
substantial. Thus the spin configuration at low temper- 
ature consists of two of the four spins in a tetrahedron 
pointing inwards, and the others outwards (the "2-in, 2- 
out" state) . Because of the statistical equivalence to the 
problem of hydrogen configuration in water ice, the 2-in, 
2-out constraint is called the " ice rule"—. Since there are 
a macroscopic number of degenerate states under the 2- 
in, 2-out ice-rule based configuration, the system exhibits 
the so-called spin ice stat o 5 : 6 : 7 with a residual entropy. 

Very recently, Katsufuji and Takagi reported the mag- 



netocapacitance of the spin-ice system Ho2Ti207^, but 
only for temperatures above 1.8 K. Here we report the 
magneto-dielectric response to temperatures low enough 
to probe the changes in the spin state associated with 
geometrical frustration. The results reveal that both the 
real part e'(T,B) and the imaginary part e"{T,B) of 
the dielectric constant can be powerful high-sensitivity 
probes for frustrated magnets. 

In the [100] field direction, in which all four spins in 
each tetrahedron are equally affected by the magnetic 
field, a crossover occurs at about 0.5 T to a state with- 
out spin degeneracy at low temperatures, in which all 
MiS point along [100] (Fig. [TJ. In the [111] field direc- 




FIG. 1: Network of tetrahedra in the pyrochlore lattice. 
Black arrows at vertices represent the magnetic moment 
rrii of Dy 3+ . Gray arrows and gray circle represent the 
total moment Mi for a tetrahedron, defined by Mi = 
X^=i m i,j- Tetrahedra (B)-(D) constitute three of the six 
nearest-neighbors for (A). 
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tion, in which one of the four spins is exactly parallel 
to the magnetic field, the 2-in, 2-out state becomes un- 
stable for fields greater than about f T and changes to 
the "1-in, 3-out state" in which all MiS have a positive 
component along the field. The pyrochlore lattice can 
be viewed as a stacking of "pyrochlore slabs" consisting 
of triangular-lattice and Kagome-lattice planes normal to 
the [111] direction. It was recently found that in an inter- 
mediate field range, the field-parallel spins on the trian- 
gular lattice align along the field, whereas the other spins 
on the Kagome lattice maintain frustration and associ- 
ated residual entropy because of the overall 2-in, 2-out 
ice-rule constraint (the Kagome-ice state) 9 ' 10 ! 11 . Because 
of such variety of frustrated spin structures which can be 
controlled by magnetic fields in particular directions, the 
spin-ice magnet is ideal to investigate how the magneto- 
dielectric response can probe the local spin structures. 



II. EXPERIMENT 

We used single crystals of Dy 2 Ti 2 07 grown with an 
infrared floating-zone furnace (NEC Machinery, model 
SC-E15HD)& and cut to a typical size of 3 x 2 x 
0.7 mm 3 with the crystalline axis [100] or [111] nor- 
mal to the wide plane of the sample. We measured 
the magneto-dielectric response as well as magnetostric- 
tion with a capacitance cell designed for a standard 
three-terminal method. We used a capacitance bridge 
(Andeen-Hagerling, model 2500A) with a fixed measure- 
ment frequency of / — 1 kHz. The capacitance cell was 
cooled down to 0.26 K using a 3 He refrigerator (Oxford 
Instruments, model Heliox VL) under magnetic fields 
of up to 7 T. The thermal contact of the sample was 
made by gold films sputtered on the top and bottom 
sample surfaces in direct contact with the copper pole 
plates, which are connected to a thermally-anchored cop- 
per coaxial cable. The gold coating was needed also to 
eliminate vacuum gap between the sample and the pole 
plates. 

The measured capacitance C was about 10 pF at room 
temperature, increased by 20 % with decreasing tempera- 
ture, and peaked at about 20 K. The measured dielectric 
loss A was of the order of 10 _1 nS at room temperature 
and decreased by up to a factor of ten at low tempera- 
tures. In all the Dy2Ti2 07 crystals we investigated, we 
reproducibly observed a strong and sharp peak in A at 
about 150 K with the width of only about 20 K. Although 
it is similar to the energy scale of the crystalline field ex- 
citation level of the Dy 3+ ion, the sharpness of the peak 
in A suggests that it is not simply due to thermal acti- 
vation of Dy 3+ moments. Its origin is not unknown at 
present. The imaginary part of the dielectric response is 
evaluated from A and C by e" = e'tan<5 = e' A/ (2irfC). 
In the data presented below, the magnetic field has been 
corrected for the demagnetization factor based on the 
magnetization data at respective temperature cfi 12 ; the 
field is expressed in the internal field B, rather than in 



the external field H. 

Since the magnetocapacitance contains information 
not only of the dielectric constant but also of changes in 
the sample thickness and area, it is important to exam- 
ine the contribution from the magnetostriction. For this 
purpose, we used the same capacitance cell but glued on 
the lower pole-plate a sample crystal coated with sput- 
tered gold film on all its surfaces. In addition, we placed 
three glass posts, each 1 x 1 x 0.84 mm 3 in size, between 
the pole plates so that there is a spacing of ca. 100 /im 
between the crystal and the upper pole plate. From the 
measurement of the capacitance with and without the 
sample crystal under magnetic fields, we extracted the 
changes in the sample shape A(d(B) / S (B)) , where d is 
the thickness and S is the area of the sample surface. The 
results for the [100] and [111] field directions are shown 
in Fig. [2] As shown below, the magnetostriction is more 
than an order of magnitude smaller than the magneto- 
dielectric response. For the [100] field direction shown in 
Fig- Eta.) , for which all the spins respond equally to the 
Zeeman energy of the external field, the magnetostriction 
increases monotonically with field. 

By contrast, we observed a clear peak for the [111] 
field direction as shown in Fig. Gib). Compared with 
the specific-heat peak position shown in Fig. |31 we con- 
clude that the magnetostriction peak reflects the change 
in the spin configuration to the "1-in, 3-out" polarized 
state. Corresponding to this change becoming a first- 
order transition below 0.4 Ki£, the magnetostriction peak 
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FIG. 2: (Color online) Magnetostriction of Dy2Ti207 for the 
(a) [100] and (b) [111] directions. d(B) is the sample thick- 
ness and S(B) is the sample area under the magnetic field B. 
Curves in (b) are offset from each other for clarity. 
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FIG. 3: (Color online) Field-temperature phase diagram of 
Dy2Ti207 for the [111] field direction. Closed triangles in- 
dicate the peak position in the dielectric loss, open squares 
the peak in the specific heafe—, and asterisks the peak in the 
magnetostriction. The lines are guides to the eye. 



is very sharp at 0.26 K. Moreover, a shoulder feature at 
0.4 T corresponds well to the change from the spin ice to 
the Kagome ice. Nevertheless, other spin-configuration 
boundaries in the phase diagram cannot be probed by 
magnetostriction within the present precision. It is in- 
teresting to note that the height of the peak remains es- 
sentially unchanged with temperature. The results here 
are semi-quantitatively consistent with those reported for 
polycrystalline Dy2Ti2C>7 at 1.7 and 4.2 K, for which the 
longitudinal magnetostriction is substantially bigger than 
the transverse oneii 

We next present the temperature dependence of the 
dielectric constant in Fig.0|and[SJ In zero field, we ob- 
served a substantial reduction in the real part e' below 
about 1.5 K, as well as a strong peak in the imaginary 



part e" at 1.2 K. These features correspond to the forma- 
tion of the local 2-in, 2-out spin configurations governed 
by ferromagnetic nearest-neighbor interaction J c ff = 1.1 
K 4 , leading to the frustrated spin-ice state with residual 
entropy. By applying a magnetic field along the [100] 
direction, the spin frustration is mostly released above 
1.0 T. Correspondingly, in a magnetic-field-cooling se- 
quence above 1 T, the features in both e' and e" observed 
at lower fields disappear completely as shown in Fig. 21 
When the field is reduced to zero after magnetic-field- 
cooling at high fields, e' smoothly decreases to the same 
value obtained in the zero-field cooling even at 0.26 K. 

We next describe the dielectric response for the [111] 
field direction, for which details of the changes in the 
spin frustration are known from the specific heat and 
magnetizationiSiii. For the [111] field direction shown in 
Fig. |S{a) , e' responds sharply to the critical end-point at 



(0.4 K, 1 T) of the first order transition from the Kagome- 
ice state with 2-in, 2-out configuration to the "1-in, 3- 
out" state. We note that the peak position of e" shown 
in Fig. EJb) shifts to higher temperature by about 0.15 
K from to 0.37 T. 

Fig. shows the field dependence of the magneto- 
dielectric response in the [111] field direction. Corre- 
sponding to the magnetization plateau in the Kagome-ice 
stated, e' displays a field-independent plateau between 
0.3 and 0.8 T at low temperatures as shown in Fig.[|Ja). 
At higher temperatures e' behaves quite differently, ex- 
hibiting a peak slightly below 1 T which broadens with 
increasing temperature. The implication of this behavior 
becomes clearer in the dielectric loss e" in Fig.[|Jb). The 
B-dependence of e" exhibits a shoulder-like anomaly at 
about 0.3 T, reflecting the change from the spin-ice to 
Kagome-ice state. The field-induced first-order transi- 
tion below 0.4 KiS results in a sharp enhancement of 
e" centered at 0.94 T. At higher temperatures e" dis- 
plays two strong peaks below and above 1 T. When plot- 
ted in the phase diagram in Fig. as closed triangles, 
these clearly correspond to the crossover lines from the 
Kagome-ice state to the "1-in, 3-random" state and to the 
"1-in, 3-out" state. In addition, another peak emerges in 
the low-field region above 1 K; the peak shifts to higher 
field with increasing temperature. This peak in e" corre- 
sponds to the crossover line between the all-random and 
the "1-in, 3-random" states. 
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FIG. 4: (Color online) Temperature dependence of (a) the real 
part e' and (b) the imaginary part e" of the dielectric constant 
of Dy2Ti207 for magnetic fields along the [100] directions. 
The data were obtained on warming after cooling under the 
presence of the magnetic field. 
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FIG. 5: (Color online) Temperature dependence of (a) the real 
part e and (b) the imaginary part e" of the dielectric constant 
of Dy2Ti2 07 for magnetic fields along the [111] directions. 
The data were obtained on warming after cooling under the 
presence of the magnetic field. 



III. DISCUSSION 

The main features in the magneto-dielectric response 
may be summarized as follows: (1) The peaks in the di- 
electric loss correspond very well to the changes in the 
spin configuration determined from the specific heat and 
magnetization. (2) The boundary between the Kagome- 
ice and "1-in, 3-random" states is precisely determined 
from e" . It is interesting to note that the boundary tem- 
perature has a maximum at 1.35 ± 0.05 K for [LqH = 0.5 
± 0.1 T. While the Kagome ice state is destabilized by 
magnetic fields for low-T, high-i? region, it is stabilized 
by fields for high-T, low-i? region. If the boundary re- 
flects the thermal activation of the least-stable spins un- 
der magnetic fields, the boundary would have a negative 
dB/dT . This is because the thermal energy to overcome 
the Zeeman energy to flip one of the spins on the Kagome 
plane with a field-antiparallcl component of the magnetic 
moment is given by 2A£q=4J c ff — (2/3)gjJ nsB and de- 
creases with increasing B. The behavior of the boundary 
at high B is explained by this mechanism. In fact, the 
above equation gives B c — 1.0 T for the disappearance 
of the first excited gap, consistent with the result shown 
in Fig. |31 On the contrary, the boundary for B < 0.4 T 
is given by Tb(B) = T si + aB with the spin-ice entering 
temperature T sl — 1.15 K and the positive slope a = 0.41 
K/T. This behavior cannot be explained even if we in- 
troduce the entropy term in the free energy and indicates 
that this boundary is essentially a collective change in the 



spin configuration, rather than by a thermal activation 
of individual spins. 

Let us now examine whether the phenomenological 
formula* relating the nearest-neighbor spin correlation 
(rrii ■ m,j) nn (here rrii = gjJi^B = 10/x_b for Dy 3+ ) to the 
dielectric constant 



e'(T,B)=e' (T)(l+a(m i -m j ) nn ) 



(1) 



is applicable to Dy 2 Ti 2 C>7. With magnetic fields in the 
[100] direction, the macroscopic degeneracy is lifted while 
the local 2-in, 2-out spin configuration is maintained. 
Thus (rrii ■ rrij} nn should be unchanged, contrary to the 
observation in which the low-temperature reduction in e' 
clearly ceases under the field as shown in Fig.^a). 

To seek an alternative phenomenological description 
applicable to Dy 2 Ti 2 07, let us replace (rrii ' mj)nn by 
(Mi ■ Mj) nn , where Mi is the total moment of the four 
spins on each tetrahedron, and the average is taken over 
the nearest-neighbor tetrahedra constituting the network 
of super-tetrahedra (shown in Fig. ^) . Such replacement 
is meaningful if a TiC<6 octahedron, which is surrounded 
by Dy 2 tetrahedra and plays the main role in the di- 
electric response 1 ^, couples more strongly with a total 
magnetization of a Dy 2 tetrahedron than with an indi- 
vidual Dy 3+ moment. Under zero field, the 2-in, 2-out 
state will yield the average (Mi ■ Mj) nn = since Mi 
randomly points in any of the six directions of (100). 
With increasing magnetic field along the [100] direction, 
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FIG. 6: (Color online) Magnetic field dependence of (a) the 
real and (b) the imaginary parts of the dielectric constant of 
Dy2Ti207 for the [111] field direction. All the data were taken 
on field reduction after zero-field-cooling and subsequent field 
application at least up to 3 T. e' has been corrected for mag- 
netostriction. See text for the dotted lines. 
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all MjS become parallel to each other and (Mi ■ Mj) nn 
= (4 x l/%/3) 2 = 16/3, if the magnitude of the original 
spin moment mj is taken as unity. In the Kagome-ice 
state under the [111] field, Mi can take only three of the 
directions and (M t • M ) nn = 16/3 x 1/3 = 16/9. With 
increasing field, the spin configuration changes to "1-in, 
3-out" state with (Mi ■ Mj) nn = 4, giving an enhance- 
ment of (Mi ■ Mj) nn by a factor of 9/4. 

Since the magnetization™ essentially saturates above 
1.5 T, we assume for such low fields that e'(T,B) is a sum 
of magnetic and non-magnetic contributions. By analogy 
to Eq. (1), we examine the relation below 1.5 T, 

e'(T, B) = 0(Mi ■ Mj) aa + A(T); (2) 

thus A(T) — e'(T, 0). In addition, for a given tempera- 
ture Ae'(T, B) = e'(T, B)-e'(T, 0) = (3(M v M ) nn . If we 
further assume that j3 is constant, Ae'(T,B) is constant 
for the Kagome-ice state since (Mi ■ Mj) nn remains con- 
stant at 16/9. Indeed, Ae'(T, B) in the Kagome-ice state 
at 0.26, 0.35, and 0.40 K shown in Fig. © (a) gives the 
same value of 0.0033, implying that (3 is about 2 x 10~ 3 . 
Moreover, in the "1-in, 3-out" state just above IT, the ob- 
served value of Ae'(T, B) is in fair agreement with the ex- 
pectation of 0.0033 x 9/4 = 0.0074 as shown in Fig.|Ua). 
Ae'(T, B) for the [100] field direction (Fig. Ufa)) is also 
in fair quantitative agreement using the same value of (3. 
Let us next consider another expression, 

a'(T,B) = /3(M*)+E'(T,0); (3) 

in which the direction u z" is defined along the field, and 
M z is the bulk magnetization : (M z = Y^li Mi Z . Mi Z 
is the z-component of Mi, N is number of spins, and 
N /4 represents the number of tetrahedra.) Incidentally, 
Eqs. 2 and 3 give the same results for the Kagome-ice 
state, the 2-in, 2-out polarized state under the [100] field, 
and the "1-in, 3-out" polarized state. However, the two 
expressions give different results for a spin-ice state under 
weak field : In the spin-ice state in the limit of low field, 
(Mi ■ Mj) nn = whereas jj(M z ) = 16/9. Furthermore, 
there is a great difference in the T-dependence. Since 
(Mi ■ Mj) nn is affected not only by the z-component Mi Z , 
but also by the transverse components Mi X and Mi y , the 
difference between (Mi ■ Mj) nn and (M z ) 2 increases with 
the increase in the thermal fluctuation, especially in zero 
and low field regions™ 

Although neither (Mi ■ Mj) nn nor (M z ) can be mea- 
sured directly, (M 2 ) can be expressed in terms of read- 
ily measurable physical quantities (susceptibility Xz an d 
mean magnetization of a tetrahedron (M z )) as (M 2 ) = 
(M z ) 2 + T\z, if thermal equilibrium is assumed. Then 
Eq. 3 can be written as 

e'(T, B) = f3{(M z f + T Xz (T, B)} + e'(T, 0). (4) 



The first term explains semiquantitatively the plateaus in 
e' in the Kagome-ice and "1-in, 3-out" states, as discssed 
above. 

The second term T\ Z {T, B) is very important to ex- 
plain the peak of e' around 1 T (Fig. HJa)), not con- 
tained in the magnetization data™. This is the region of 
metamagnetic instability with large Xz- The factor T in 
the second term should make the contribution smaller at 
lower temperature, while the factor Xz makes the peak 
broader at higher temperatures, just as observed. Equa- 
tion 4 thus captures the information of e' which contains 
not only the static magnetization (M z ) but also the mag- 
netic susceptivility Xz- 

A preliminary Monte-Carlo simulation suggests that 
the observed peak in e' at about 1 T is much greater than 
that expected from eq. 4 and the expression in terms of 
(Mi ■ Mj) nn better describes the observation. Yet addi- 
tional contribution seems to be present at higher fields: 
even for the case where the magnetization is saturated 
with increasing field, the dielectric constant keeps chang- 
ing. The sign of the change is negative in the [100] field 
direction above 2 T (Fig.JSfa)), while it is positive in the 
[111] field direction. This difference in sign may provide 
an important clue to understanding the variation of e' at 
high fields. 



IV. CONCLUSION 

We have demonstrated that magneto-dielectric re- 
sponse can be used as a high-sensitivity and high- 
precision probe with fast data acquisition capabilities for 
the investigation of a geometrically frustrated magnet. 
We found that the dielectric loss e" is particularly sen- 
sitive to probe subtle changes in the spin configuration, 
not easily detectable by magnetostriction. In addition, 
the dielectric constant s' contains information of not only 
the magnetization (M z ), but also of the susceptibility Xz- 



Acknowledgments 

The authors are grateful to T. Katsufuji for useful 
discussion and sending their data prior to publication, 
to Koichiro Tanaka for letting us use a sputter coating 
machine, to M. Gingras for useful discussion and infor- 
mation, and to R. Perry, K. Deguchi, K. Kitagawa, K. 
Ishida, H. Yaguchi, S. Nakatsuji, and S. Grigera for their 
technical help as well as useful advises. This work is sup- 
ported by Grants-in-Aid for Scientific Research (S) from 
JSPS and for the 21 Century COE program "Center for 
Diversity and Universality in Physics" from MEXT of 
Japan. 



6 



* Electronic address: msaito@scphys.kyoto-u.ac.jp 

1 for example: A. P. Ramirez, Handbook of Magnetic Mate- 
rials 13 (ed. K.H.J. Buschow, Elsevier Science B.V., 2001) 
Chap. 4. 

2 T. Katsufuji and H. Takagi, Phys. Rev. B 64, 054415 
(2001). 

3 T. Kimura, T. Goto, H. Shintani, K, Ishizaka, T. Arima, 
and Y. Tokura, Nature 426, 55 (2003). 

4 S.T. Bramwell, and M.J. P. Gingras, Science 294, 1495 
(2001). 

5 A.P. Ramirez, A. Hayashi, R.J. Cava, R. Siddahthan, and 
B.S. Shastry, Nature (London) 399, 333 (1999). 

6 H. Fukazawa, R. G. Melko, R. Higashinaka, Y. Maeno, and 
M.J. P. Gingras, Phys. Rev. B 65, 054410 (2002). 

7 J. Snyder, B.G. Ueland, J.S. Slusky, H. Karunadasa, R.J. 
Cava, and P. Schiffer, Phys. Rev. B 69, 064414 (2004). 

8 T. Katsufuji and H. Takagi, Phys. Rev. B 69, 064422 
(2004). 

9 K. Matsuhira, Z. Hiroi, T. Tayama, S. Takagi, and T. 
Sakakibara, J. Phys.: Condens. Matter 14, L559 (2002). 



Z. Hiroi, K. Matsuhira, S. Takagi, T. Tayama, and T. 
Sakakibara, J. Phys. Soc. Jpn. 72, 411 (2003). 
R. Higashinaka, H. Fukazawa, and Y. Maeno, Phys. Rev. 
B 68, 014415 (2003). 

T. Sakakibara, T. Tayama, Z. Hiroi, K. Matsuhira, and S. 
Takagi, Phys. Rev. Lett. 90, 207205 (2003). 
L. G. Mamsurova, K. K. Pukhov, N. G. Trusevich, and L. 
G. Shcherbakava, Sov. Phys. Solid State 27, 1214 (1985). 
R. Higashinaka, H. Fukazawa, K. Deguchi, and Y. Maeno, 
J. Phys. Soc. Jpn. 73, 2845 (2004). 

M. Itoh, R. Wang, Y. Inaguma, T. Yamaguchi, Y-J. Shan, 
and T. Nakamura, Phys. Rev. Lett. 82, 3540 (1999). 
Instead of using the bulk magnetization (M 2 ), one may in- 
troduce thermal average of the square of the local magne- 
tization along the applied field (M 2 ). However, this quan- 
tity neither gives zero for the low- field limit, nor can be 
expressed in terms of other measurable bulk quantities. 
Thus we use (M 2 ) in the following discussion. 



